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Abstract 

We study the convergence rate of the solutions of the incompressible Euler-a, an inviscid second-grade 
complex fluid, equations to the corresponding solutions of the Euler equations, as the regularization pa- 
rameter a approaches zero. First we show the convergence in H s , s > n/2 + 1, in the whole space, 
and that the smooth Euler-a solutions exist at least as long as the corresponding solution of the Eu- 
ler equations. Next we estimate the convergence rate for two-dimensional vortex patch with smooth 
boundaries. 



On 

1 Introduction 



The equations of motion for a visco-elastic second-grade non-Newtonian complex fluid are given by the 
system (see, e.g., [27,28,72]) 



i/Au+(u-V)v + VjVuj + Vp = 0, (1.1) 



dv 

v = (l - a 2 A) u, 
V-u = 0, 

v(x,0) =v in (x), 

where the fluid velocity held, v, and the pressure, p, are the unknowns; v m is the given initial velocity held, v 
is kinematic viscosity, a > is a material parameter which represents the elastic response of the fluid, and we 
use Einstein's summation convention. The inviscid version of this model, i.e., when v — 0, is mathematically 
identical to the Euler-a [also known as the Lagrangian-averaged Euler-a] model, which was independently 
introduced and derived in the Eulcr-Poincare variational framework in [36,37]. In this variational theory the 
parameter a is interpreted as a spatial filtering scale of the velocity field v. 
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In [17-19, 29, 30] the corresponding Navier-Stokes-a (NS-a) [also known as the viscous Camassa-Holm 
equations or the Lagrangian-averaged Navier-Stokes-a (LANS-a)] model, which is a regularization of the 
Navier-Stokes equations (NSE), was obtained by introducing an appropriate ad hoc viscous term into the 
Euler-a equations, that is, by adding the viscous term — vAv, instead of — i/Au in (JTTTJ) . While the question 
of global regularity of three-dimensional (3D) viscoelastic model (|1.1|) is still a challenging open problem, 
the 3D NS-a model is proven to be globally well posed [30]. The extensive research of the a-models (see, 
e.g., [2, 4,5, 9, 12-14, 17-23,29-31, 31-33, 38,39,43,44,49,50,52, 55,62, 65, 74]) stems, on the one hand, from the 
close agreement of their steady state solutions to averaged empirical data, for a large range of huge Reynolds 
numbers, for turbulent flows in infinite channels and pipes [17-19]. On the other hand, the a-models, for small 
values of the parameter a, can also be viewed as numerical regularizations of the original, Euler or Navier- 
Stokes, systems. The main practical question arising is that of the applicability of these regularizations to 
the correct predictions of the underlying flow phenomena. In particular, it becomes important to investigate 
the problem of convergence of the a-models, as the regularization parameter a approaches zero. This 
problem has been studied in various contexts. In [30] the convergence, as a — * 0, of a subsequence of the 
weak solutions of the three dimensional (3D) NS-a equations to a Leray-Hopf weak solution of the 3D NSE 
equations is shown, for the case of periodic boundary conditions. Similar results are also reported in [55] 
concerning the MHD-a model. In [21,75] it is shown that the trajectory attractors of the 3D Leray-a and 
NS-a models converge to the trajectory attractor of weak solutions of the NSE as a — > 0. The authors 
of [16] obtained a rate of convergence of the solutions of the 3D NS-a equations with periodic boundary to 
the solutions of the Navier-Stokes equations, as a — > 0, for small initial data in Besov-type function spaces, 
for which global existence and uniqueness of solutions can be established. Recently, the convergence rates of 
solutions of various two-dimensional (2D) a-regularization models, subject to periodic boundary conditions, 
toward solutions of the exact Navier-Stokes equations, as a — > 0, have been studied in [14]. For 2D Euler-a 
regularization the authors of [5] show the convergence of a subsequence of the weak solutions of the Euler-a 
equations with a distinguished sign vortex sheet initial data (a Radon measure supported on a curve) to a 
solution of the 2D Eulcr equations, as a — ► 0. We elaborated on this result below. It is also worth mentioning 
that the problem of weak convergence of solutions of the viscous second-grade equations (jl.ip with L 2 initial 
data to a solution of the Navier-Stokes equations, as a — > 0, is treated in [41,42]. 

In this paper we follow the above mentioned philosophy proposed in [55] and consider the Euler-a model 
as a numerical inviscid regularization of the Euler equations. We study the convergence rate of the solutions 
of the 2D and 3D Euler-a equations to the corresponding solutions of the Euler equations for smooth initial 
data in the whole space, as the regularization parameter a approaches zero. In the 2D case we also investigate 
the convergence rate of the solutions of the Euler-a to the corresponding solutions of the Euler equations 
for vortex patch initial data. This program is analogues to the established results concerning the rate of 
convergence of the NSE to the Euler equations, as the viscosity v — > 0, see [1,26,45,60]. However, in the 
2D case there is an advantage of the Euler-a regularization over the Navier-Stokes equations regularization 
since the former regularizes the solution by transporting the vorticity with a smooth vector field, and hence 
it preserves the structure of the vortex patch and vortex sheet while regularizing the motion. 

The incompressible Euler equations are given by 



where v, the fluid velocity field, and p, the pressure, are the unknowns, and v m is the given initial velocity 
field. For results concerning the Euler equations, see [57,58], and for recent surveys, see [6,25]. 

In M 2 , the vorticity formulation of Euler equations is obtained by taking a curl of (II. 2| and is given by 



— + (v ■ V)v + Vp = 0, 
at 

V • v = 0, 

v(x,0) =v m (x), 



(1.2) 




(1.3) 
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q(x,0) = q in (x), 

where K (x) = ^rV 1 " log |x|, q — curlu is the vorticity, and q m is the given initial vorticity. 

Yudovich [78] obtained the existence and uniqueness of weak solutions of the 2D incompressible Euler 
equations for initially bounded vorticity (see, also, [3, 48] for an alternative proof, [73] for an improvement 
with vorticity in a class slightly larger than L°°, and [68] for review of relevant two-dimensional results). In 
particular, the problem of evolution of vortex patches, where the vorticity is a multiple of the characteristic 
function of a bounded domain, has a unique global solution, and it was proved in [15] (see also [10]) that 
C 1 ' 7 , 7 > 0, boundaries of the patches remain C 1 ' 7 for all times. In [26] it was shown that the L 2 norm 
of the difference between the solutions of NSE and the corresponding solution of the Euler system for such 
initial data converges to zero, as the kinematic viscosity v — * 0; even though none of the solutions is in L 2 . 
The (vt) 1 ^ 2 rate of convergence of [26] was improved to (vt) 3 ^ 4 in [1], due to the fact that the vorticity of 
the vortex patch with C 1 ' 7 boundary is in fact in a Besov space -B^^, see also [60] for a simpler proof and 
an extension to K 3 . 

As we have mentioned above the Euler-a model [19,24,35-37,59] is an inviscid regularization of the Euler 
equations (|1.2[) . which is given by the system 

dv a 

— + (u a - V) u Q + u"V< + Vp a = 0, (1.4) 
at 

v a = (l-a 2 A) u a , 
V • u a = V • v a = 0, 
v a (x,0) = v in ' a (x), 

Here u a represents the "filtered" fluid velocity vector, p a is the "filtered" pressure, a > is a regularization 
length scale parameter representing the width of the filter. Observe that for a = one recovers, formally, 
the Euler equations (|1.2|) . The vorticity of the 2D Euler-a model q a = curlw" obeys the equations 

?£ + (u a -V)q a = 0, (1.5) 
/ =K a *q a , 
q a (x,0) = q in > a (x). 

The smoothed kernel is K a — G a * K, where G a (x) — -^z-^Kq (J^j is the Green function associated 

with the Helmholtz operator (/ — a 2 A) (see, e.g., [66]), the function Kq is a modified Bessel function of the 
second kind of order zero (see, e.g., [77]). 

The 2D Euler-a equations were studied in [63] , where it has been shown that there exists a unique global 
weak solution to the Euler-a equations with initial vorticity in the space of Radon measures on R 2 , with a 
unique Lagrangian flow map describing the evolution of particles. We remark, however, that the question of 
global existence of weak solutions for the three-dimensional (3D) Euler-a equations is still an open problem. 
In [44] the global existence of weak solutions is shown for the 3D axisymmetric Euler-a equations without 
swirl, for initial vorticity being a finite Radon measure with compact support, also, the global existence and 
uniqueness is established for the compactly supported vorticity in L p , p > 3/2, see also [11]. 

As in the 2D Euler case, the vortex patch is transported by the flow under evolution of Euler-a equations. 
The following result is the essence of Theorem IA.ll stated below. If the initial vorticity is a multiple of the 
characteristic function of a simply connected bounded domain f2 m , with a certain technical condition on the 
boundary being a simple curve, and the boundary 9f2 m is in either one of the following spaces: Lip, or C 1 '^ , 
< (3 < 1, or C 2 '' 3 , < (3 < 1, or C n '^ , n > 3, < /3 < 1, then the boundary of the vortex patch remains in 
the corresponding space for all time. 

The outline of the paper is as follows. In Section [2] we study the convergence rate of the solutions of 
the Euler-a equations to the solution of the Euler equations for strong solutions that belong to the Sobolev 
space H m (K™), m > n/2 + 1, for n = 2,3. We show that interval of existence of Euler-a solutions contains 
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the interval of existence of the corresponding Eulcr solution, and that the convergence is uniform in time, 
for time intervals compactly contained in [0, T*), where T* is the time of existence of the solution of the 
Euler equations. We also show that in the H m ~ 2 norm the solution of the Euler-a equations differs from the 
solution of the Euler equations by order a 2 . It is worth mentioning that this result corresponds to the vt 
convergence rate of the solutions of the NSE to the one of the Euler equations, for the inviscid limit of the 
classical solutions of the NSE equations in the whole space see [45,46,60,61,67]. The issue of inviscid limit 
of the NSE in domains with physical boundaries, subject to the no-slip Dirichlet boundary conditions, is a 
very important open problem, for both theoretical study and applications. The problem emerges first from 
the boundary layer, which appears because we can not impose a Dirichlet boundary condition for the Euler 
equation, then the nonlinear advection term of the Navier-Stokes equations may propagate this instability 
inside the domain. Very few mathematical results are available for this very unstable situation. One of the 
most striking results in this direction is a theorem of Kato [47], see also [69,76]. 

In Section[3]we study the convergence for the vortex patch problem. Specifically, we show the convergence, 
as a — » 0, of the L 2 norm of the difference between the solutions of the Euler-a equations and the solution 
of the Euler equations for the vortex patch initial data with vorticity being a characteristic function of a 
simply connected bounded domain with C 1,7 , 7 > 0, boundary, even though neither of the solutions are 
in L 2 . The convergence rate is of order (a 2 ) 3 ^ 4 , which corresponds to the optimal convergence rate of the 

difference between the solutions of the NSE and the Euler equations which is of order (vt) 3 ^ 4 [1]. 

We remark that, ideally, we would like to compute the rate of convergence in dimensionless units, however 
due to the absence of typical length scale in R n , the above rates of convergence of the NSE to the Euler 
equations are given as powers of (vt), which has the units of length square, see, e.g., [1,26,45,60]. Similarly, 
in our case, the rate of convergence involves the parameter a 2 . Observe, that one can artificially cook up a 

length scale from the initial value , e.g., ||v m || i2 ||9 OT | L2 i n the case of H m (R™), m > n/2 + 1, solutions, or 
11 1 1 1/2 11 1 1 — 1/2 

l|9 m |lii In m ^ ne 2D vortex patch case. Contrary to this, if one is interested in bounded domains or 

domains with periodic boundary conditions, typical length scale will be dictated by the size of the domain, 
and then the rated of convergence will be expressed as dimensionless quantities. 

In the following all the constants C are independent of a, and all the a dependencies are spelled explicitly. 

2 Classical solutions 

In this section we study the convergence rate of the solutions of the Euler-a equations to the solution of 
the Euler equations for strong solutions that belong to the Sobolev space H m (R n ), m > n/2 + 1. For 2D 
incompressible flow with initial velocity v m S H m (R 2 ), m > 3, the unique solutions of the Euler equations 
exist globally in time (see, e.g., [57]). Similarly to the Euler case, the Euler-a equations has a unique global 
solution, since, as in the 2D Euler case, we have an a priori uniform control over the L°° norm of the 
vorticity, which implies the global existence, as in the proof of the Beale-Kato-Majda criterion [7]. In M 3 
only local in time existence of strong solutions of the Euler equations has been shown, see [45,46,53,54,57]. 
The existence and well-posedness of the Euler-a equations for a short time can be easily shown following the 
classical theory of the Euler equations; see also [40, 56] for an analogue of the Beale-Kato-Majda criterion 
for the Euler-a model. More precisely, one has the following result: 

Proposition 2.1. Let v m e H m (R n ), m > n/2+1. There exists T* = T* (\\v m \\ Hm ), T* > ^Jjj - , such 

that for any T < T* there exists a unique solution v £ C([0,T];H m (R n )) n AC ([0, T]; iJ m_1 (R™)) of the 
Euler equations (|1.2[) with initial data v m . In two dimensions the solution exists globally in time. Similar 
results hold for the Euler-a equations (jl.4p with the maximal interval of existence of the three-dimensional 
Euler-a equations being also dependent on a. 

In the next theorem we show that the solutions of the Euler-a equations for the H m , m > n/2 + 1, initial 
data, exist at least as long as the solution of the Euler system exists, and converge, as a — > 0, to the solution 
of Euler equations. Our proof of this result follows the ideas in [60]. 



4 



Theorem 2.2. Let v in ,v in > a £ H m (R n ), m > n/2 + 1, and \\v in - v in ' a \\ Hm -> 0, as a -» 0. Lei T* &e i/ie 
time of existence of the solution of the Euler system (TO]) v £ C; oc ([0, T*);H m (R n ))nAC; oc ([0, T*); iJ m_1 (IT)) 
wzi/i initial data v m . Then, for all < T < T* f/iere exists < a = a (v m , ir* 11 '", J 1 ) suc/i i/ia£ /or all a < a 
there is a unique solution v a £ C([0, T];H m (W 1 )) n AC ([0, T];i/" m_1 (M n )) o/ tfie LuZer-a equations (fl~4j) 
tuif/i initial data v m ' a . Moreover, 

0, 



lL°°([0,T],ff m ) 



as a — > 0, and /or all <t <T 

\\(v a -v) (*)H« — » < (|K 

(t)llfl»-l < (IK 

where C = C (||^ m || ffm j L) *s independent of a. 
We use the standard Lf m (R n ) norm denned by 



1^ 



Iff" 



Ca 2 i) e ct , 
Cat) e ct , 



(2.1) 
(2.2) 



where / denotes the Fourier transform of /. 

To prove the theorem we need the following estimates (see, e.g., [45,58]) 

Lemma 2.3. Let n = 2, 3. Let u, v £ H m (R™) , divu = 0, then there exist a constant C > 0, depending on 
m, such that 



\((u ■ V)v,v) Hm | < C \\Vu\\ Hm ^ \\v\\ 2 Hm < C \\u\\ Hn \\vf Hm , m > f + 1, 
|((u ■ V)v, v) Hm | < C \\Vu\\ H2 \\v\\ 2 Hm < C \\u\\ H3 \\vf Hm , 



(2.3) 



Let u £ H m (R n ), v £ H m+1 (R n ), and Zei ^ (u, u) 6e either one of the following bilinear forms: \& (u, u) = 
it x (V x v), ^ (u,v) — (u ■ V)w or ^ (u,w) = v.Wiij, then there exists a constant C > 0, depending on m, 



such that 



and for m > \ + 1 



\\*(u,v)\\ Hm <C\\u\\ Hm \\v\\ H 

\\*(u,v)\\ Hm <c\\u\\ Hm \\v\\ H3 , 



m > 



+1 j "« 2 : 



||*(«,«)||^<C(||u|| ffm ||«|| Jm + ||u|| io 
We also use the following lemma 
Lemma 2.4. Let g £ L 2 (R n ) and / = (l - a 2 A) _1 .g. Then 

(-A^Vl^-llfflli*. 



NU+0- 



(2.4) 
(2.5) 



(2.6) 



Proof. By taking a Fourier transform we have 

/(£) = 



3(0 



i + « 2 i^r 



hence 



(-A) 1/2 / 



L 2 (R") 



19 (01 



(i+« 2 icf 



< 



||3|| i 2 (K „ ) sup 



y>0 (1 + y) 



1 

^4 



□ 
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Next, we prove Theorem l2.2l 

Proof. First, let us assume that v,v a e C([0,T);H m (W))nAC ([0, T]; H m - 1 (M™)) are the solutions of the 
Euler and the Euler-a systems (ll.2p and (|1.4I) with initial data w m and v m ' a , respectively, on some mutual 
time interval [0, T], with 

IML-acT],^) , ll« a IL=o a o,T],^) < c (r, 1^1^) , (2.7) 

and let us show the convergence rates in H m ~ 2 and H m ~ l . In the second part we show that the solutions of 
the Euler-a equations exist at least up to the time of existence of the solution of Euler equations and satisfy 

(E3). 

The difference w a — v — v a satisfies the following equation 

^ + {w a • V)« + VjVwf + (u a • V)w a + w Q V«? (2.8) 

- (a 2 Au a ■ V)v - VjV (a 2 Au^) - VjVvj + V(p-p a ) = 0, 
w a (x,0) = v m {x,0) - v m - a (x,0). 

Let k be either m — 1 or m — 2. Taking the H k inner product of (|2.8[) with tc" (t), we obtain 



where 



Ji = \(w a x (V x u),u; Q ) Hfe | 
J 2 = |(( M «-V) W Q , W Q ) fffc |, 

/ 4 = Q 2 |(A/x(Vx v ),w a ) flt |, 



due to the identity 
Now, by |23J) 
by (E3]) 
by (1231) 
and by (|2~4]) 



(6 • V) a + ajVbj = -b x (V x a) + V (a • 6) . (2.9) 



/i<c|M| Hm |K||^, 
/ 2 <c|KIU. IKII 2 H ,<cK|| ffm |K||^, 

h<C\\u^ Hm \\w^ 2 Hk <C\\v a \\ Hm \\w«f Hk , 

I 4 <a 2 \\Au a \\ Hk \\v\\ Hm \\w a \\ Hk . 
Summing up we get the energy estimate 

~ \\w a f H * < C \\w a f Hk (\\v\\ Hm + \\v a \\ Hm ) + a 2 \\Au a \\ Hk \\v\\ Hm \\w a \\ Hk . 
Now, for k — m — 2 we have that 

\\Au a \\ Hm _ 2 <C\\u a \\ Hm <C\\v a \\ Hm , 

and hence 

™ \\w a \\% m - a < C \\w a f Hm _ 2 (\\v\\ Hm + \\v a \\ Hm ) + Ca 2 \\v a \\ Hm \\v\\ Hm \\w a \\ Hm _ 2 
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while for k = tyl — 1 we have by (|2.6p 

a\\\u a \\ Hm ^ <C||(-A) 1/2 « Q 

and hence 



„ ^ciKIU, 



^ ll«^llir™-x < C II^H^-x (||«||^ + H^ll^™) + C« ll^ll^™ NIh-KV 



Therefore, by Gronwall lemma and (|2.7|) we obtain that 

IK < (IK (• > 0)|| flW _ a +Co a *) e ct , 

and 

< (IK 0, 0)11^ + Cat) e ct , 

where C = C (\\v m \\ Hm ,T). 

Next, we show that the solution of the Euler-a equations (| 1 .4|) exists as long as we can solve the 
Euler system (|1.2|) . and that v a converges to v in L°° ([0, T] , if™ 1 ), as a — > 0. The proof follows the 
ideas in [60]. We regularize the initial data by taking v m > s = J 7-1 (x|f|<i/<5 (0 F (v m )) and v m > S ' a = 
F^ 1 {X\i\<i/S {£)F{v in,a )), for some 8 e (0,£ ]. Let a* be such that v m > a " ^ (if v in > a = for all 
a, then the proof is trivial). Since ||w m,a — v m \\ Hm — > 0, as a — > 0, then there exists ao such that 
K n ' a I rrr „ < u* n | ,,_ + ||w m,a II „_ for all a < an. The regularized initial velocities of the Euler and 

II \\ rl II 1 1 tt II Mil 

of the Euler-a equations satisfy, for all a < ao, 

II in, a, 5 II \\ v in - S \\ < A" 

lr \\H m ' II llff" - ' 

II ■ All ,i ■ All K 



and for s£ [0, ; 



v m ' Q ' 5 || ffs , ||w m - v m > 6 \\ Hs < K5 rn ' s , 



where K = K (||u m ||„ m , ||w m ' a || ffm )- In the following we fix s such that § < s < m— 1. 

Let i>, i/, u", v a ' S be the corresponding solutions of Euler and Euler-a equations with initial data v m 
and v m ' S , respectively. Notice that, as explained below, the solutions v, V s , w Q and /■* exist on some time 
interval [0, To], Q < Tq < T* , which is independent of a and S, and also, for all t £ [0, To] 

\\v(t)\\ H ™,h 5 (t)\\ Hm <C(T ,K), \\<P S (t)\\ Hm+1 < C{To s ' K) , (2.10) 

where y denotes v or u Q , and tp s denotes either v s or v a ' S . Indeed, in R 3 , writing Euler and Euler-a equations 
in the vorticity formulation we have that q — curl v and q a = curl v a satisfy 

^ + (« • V) q = {q • V) V 

and 

respectively. Making the H k estimates (one can use (|2.3j) and the fact that for k > 3/2, H k (R 3 ) is a Banach 

3 
2 : 



algebra), due to m — 1 > I, we obtain 



^IIV-IU-! <c|HI 2 ffm -x, 

All,/-- 5 !! ^ ^ IL/.5I' 2 

eft 
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where ip denotes q or q a , and ip s denotes q 5 or q a ' S , and hence the solutions exist on a certain interval [0, To] 
depending only on K and 8q, independent of a and 6. Furthermore, we haveQ 



•(t)\\ H ^<C(T ,K),\\^(t)\\ Hm < 



C(T ,K) 



Hence, from the Biot-Savart law, ip — J R3 |^_Jja x curl (,5 (y) dy, we obtain (|2.10p . In R 2 the solutions of 
Euler and Euler-a equations exist for all times, see Proposition ^. II 

Now, one can show, following [60], that v s is a Cauchy sequence in the Banach space C([0, Xb] ,H m ), 
due to 



v s - v 5 ' 



< 



v in,5' _ y in,S 



+ (max{ ( 5, ( 5'}) m - s - 1 C(ro,|h m || ffm ))e C ( To 'll"" 1 |l^). 



C([0,T ],ff"») 

The limit v of v s ,as 8 — > 0, is the solution of Euler equations, and from (|2. 12|) . we have 



(2.12) 



\v — v 



lL°°([0,T( 



h]Bm) < (II^-^'IIh* + C (T ,\\tf n \\ Hm ) S™- 8 ' 1 ) Ce c ( To '\\ vin \\*™\ (2.13) 



Next we show that v a - s is a Cauchy sequence in C ([0, To] , H m ), and converges to the solution of Euler-a 
equations v a , as 6 — > 0, and we have 



\v a -v a ' 5 \ 



L°°([0,T< 



•o],*») " dl U 



\ Hm +C(T ,K)6 m - s - 1 )Ce c( - To - K l (2.14) 



Specifically, in R 3 (R 2 case can be done similarly), we first assume that S > 6', and denote zu a ' S ' 5 

q a,S' _ q a,S^ w a,S,5' _ ^,5' _ ^o,^ We have that 



dvu 



a.S.S' 



-{u a ' 5 ' 



V) zn a > 5 ' S '+( ( (1 - a 2 A) L w a > 5 ' 5 ' ) • V) q 



dt 

Using (12. 3|) and ()2.5|) we obtain 

2 



( 'a,S,S' 
[ZD 



W)u a ' 5 ' + (q a ' s -V) ((l-a 2 A)" l tf 



< c 



a, 6' 



Iff""! 



a, 5, 5 



-c 



la Q,,5 ll 
|y I Ih" 



a, 5, 5' 



Now, the difference w a ' S ' S = v a,s — v a ' S satisfies the equation (see (|2.9[) 



dt 



+ ((1 - a 2 A) 1 x (V x iA 5 ') + x (v x io°' y ') + V (p a ' S ' - p a ' S ^j = 0, 



where p a , p a ' S are the modified pressure. By (|2.3|) and (12. 4p we obtain that 

< C* ( \\v a - A \\ TI + v a - d ' 

H 



1 d 

2dt 



:,S,S' 



IN>(*)llffm-l <*(«) 



J It follows that 

(2.11) 

where \E> is the solution of the scalar value initial-value problem (t) = C^ 2 (t), where (0) is the H m ~ l norm of either 
one of q in , 5*"'°, <j in -' 5 , q" 1 ^.*, * (0) < X. 

Now (t) = 1 _gt^(Q) exists on a certain interval of time [0, To], To > Q^fm > cTk' ODV i° us ly independent of a and 8, and 
we have sup 0<t<T() \& (t) < t _ ^ ^ . Now, if a solution ip (t) existing on [0, T Q ' 4 ] , such that T a ' S < To, then the system can be 
solved with initial value ip 5 (T a - S ) £ H m ~ x , to continue the solution to [0, T a - S + T°'*l in which (f2TTTb is true. Iterating this 
argument we can continue the solution to cover the whole interval [0,T] with the estimate l|2.11|l throughout. (If the solution 
cannot be continued at some time f < To, then limsup t _ ( y_ \\ip s {t)\\ H m-i = °°; a contradiction to 112.111 1.') 
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Recall, that j < s, hence, by Sobolev embedding theorem and Gronwall lemma we have 



AS' 



L°°([0,T o ];L° 



< 



< 



,,Oi,5,8 f 



w a > s ' 5 ' (0) 



L=°([0,T};H") 

e c ssb a, '^)\\ am +\\^-''(T)\\ Km dr 



< C(T a ,K)5 m - s e c{T «- K l 



Using also that ||<J Q ' ,5 || Hm < C(T ° 5 ' K) we obtain 



d_ 

dt 



zu 



:A5' 



<G\ 



H" 



H" 



+ Q 



w ■ 



AS' 



H" 



+ C(T ,K) g*->-i e cm>,K) 



hence 



I 



lL°°([0,T, 



b];H— i) - ( 



Hmi +S m - s - 1 C(T Q ,K))e 



C(T ,K) 



Making the same estimates for the case 8' > S, we obtain that 



~ a ?s I 



lL°°([0,T ( 



iliff™- 1 ) " ( 



Hmi +(m S K{S,S , }) m - s - 1 C(T ,K))e c ^ K \ 



which implies that w a,s is a Cauchy sequence in the Banach space C ([0, To] , H m ), and its limit v a , which 
is the solution of Euler-a equations satisfies (|2 . 14[) . 

It remains to show that \\v s (t) — v a ' S (t) „ m converge to zero, as a and 5 converge to zero uniformly in 
[0, T ], and then the convergence of v a to v in L°° ([0, T ] , H m ) follows. 

From (|2.8[) and Lemma \2. 31 we have that the difference w s — V s ~ v a - 5 satisfy 



^ " ,.<5|l ^ /IL.^II i IL, a :^ll \ IL.,^11 i /~i IL.,^11 /IL.^II i IL,c*,<5| 



dt 



\w s \\ Hm < C (\\v s \\ Hm + \\w s \\ Hm + C \\w s \\ Lao (\\v s \\ Hm+1 + 

+ Ca 2 \\Au a ' S \\ Hm \\v 5 \\ Hm + Ca 2 \\ Au a ' S \\ Lae \\v s \\ Hm+1 . 



\H m + 1 



By Sobolev embedding theorem, (|2.2p and Gronwall lemma we have that 



IZ,o°([0,To],L°°(R")) 



-i, < {\\v m ' S - u in,Q ' 4 || ffm _ 1 +C(T 0l K)a) e c{To ^ 



li°»([0,T Q ],H w '- 1 ) 

< (\\v m -v m ' a \\ Hm ^ +C(T 0l K)a) e c ( ;:! k ' 



by 



Hence 



a\\Au a \\ Hm <C\\v a \\ Hm+1 , 



<a\\Au a - s \\ Hm ^<C\\v a \\ H , 



-\\w s \\ Hm <C\\w s \\ Hm +C 



e° + C- 



S 5 18 

where C = C (K, To). By using Gronwall lemma, first letting a — > 0, and then letting 8 — > 0, while choosing 

0, we obtain that \\v s - v a < 5 \\ raofln ^ nm ^ -> 0. Using also (j2~T3)) and ([2l4]) 



5 such that 



s ■ ~, — — — !r - lli°°([o,r ],H m ) 

the convergence of v a to u in L°° ([0, To] , iJ m ) follows. 

Now, let the Euler solution exist on [0, T*), then we can continue the solution of the Euler-a equations up 
to any T < T* in a finite number of iterations using the above argument. Indeed, we continue the solution 
from a time interval [0, To + ... + Tt_i] to a time interval [0, To + . . . + Tc], by solving the Euler-a equation 
with the initial data v a (To + . . . + Tt_i), a < min{ao, . . . , ak-i}, which converges to v (To + . . . + Tk-i) 



for a time T k > C (\\v (T k ^)\\ Hm + \\v m > a * \\ Hm ) 1 > C (\\v\ 



L°°([0,T},H m ) 



□ 
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3 The vortex patch case 



In this section we study the convergence, as a — > 0, of the L 2 norm of the difference between the solutions 
of the 2D Euler-a equations (|1.4[) and the solution of the 2D Euler equations (|1.2p for the vortex patch 
initial data with vorticity being a characteristic function of a simply connected bounded domain with C 1 ' 7 , 
7 G (0,1), boundary. We show that the convergence rate is of order (a 2 ) 3 ^ 4 , which corresponds to the 
optimal convergence rate of the difference between the solutions of the 2D NSE and the 2D Euler equation 
which is of order (vt) 3 ^ 4 . 

Yudovich [78] obtained the existence and uniqueness of weak solutions of Euler equations for initial 
vorticity in L°° (R 2 ) n L 1 (R 2 ), in particular, for the problem of evolution of vortex patches, where the 
initial vorticity q m is assumed to be proportional to the characteristic function of a bounded domain f2 m , 
q m = qoXn in - Due to the conservation of the vorticity along particle trajectories, the vorticity q (t) remains 
a characteristic function of an evolving in time domain f2 (t) . For the case where the boundary of the patch 
n m belongs to C 1 ' 7 , 7 > 0, it was proved in [15] (see also [10]) that the Euler system has a unique solution 
v 6 Lf£ c (R, Lip (R 2 )), and f2 (t) remains a bounded C 1,7 domain. It was proved in [26] that the L 2 norm 
of the difference between the solutions of NSE and the corresponding solution of the Euler system for such 
initial data converges to zero, as the kinematic viscosity v — > (even though none of the solutions are in 
L 2 ). The rate of convergence was improved to (vt) 3 ^ 4 in [1], due to the fact that the vorticity of the vortex 
patch with C 1,7 boundary is in fact in a Besov space B^^, see also [60] for a simpler proof and an extension 
to R 3 . 

It has been shown in [63] that there exists a unique global weak solution to the Euler-a equations for 
initial vorticity in yV((R 2 ), the space of finite Radon measures on R 2 , with a unique Lagrangian flow map 
describing the evolution of particles. In [5] we show the convergence, as a — > 0, of the weak solutions of 
Euler-a equations with a distinguished sign initial vorticity in .M(R 2 ) n H^, (R 2 ) (vortex sheet data) to 
those of the 2D Euler equations. 

Since the solution q a of the Euler-a equation is transported by the smoothed vector field u a , then in 
the Euler-a case the vortex patch is also transported by the flow. The global existence and uniqueness 
results for the smooth vortex patch evolution under the Euler-a equation can be obtained using arguments 
similar to those presented in [5] and [57, Chapter 8]. Specifically, if q m (x) — qoXn in i x ) 1S a multiple of the 
characteristic function of a simply connected bounded domain Q ln , and the boundary dft m is in either one 
of the following spaces: Lip, or C 1 ' 13 , < (3 < 1, or C 2 ' 13 , < j3 < 1, or C n < , n > 3, < < 1, then the 
boundary of the vortex patch remains in the same space as dVl m for all times. We describe this result in 
details in Appendix, Section lA.il 

In studying the convergence rate we use the following results. 

Proposition 3.1. Let q ln 6 L°° (R 2 ) n L 1 (M 2 ) . Then there exist unique global solutions q and q a of Euler 
and Euler-a equations (|1.3p and (|1.5p respectively. Moreover, the L p norms of q (of f)l .3(1 ) and q a ( of (|1.5p ) 

are conserved, namely, \\q (•, t)\\ LP = \\q a (•, i) \\ LP — ||9 m || iP jl < P < °°- I 71 addition, the velocities are 
bounded uniformly 

\\v(;t)\\ L ~,\\v a (;t)\\ L ~ <(ll<rlUk m ILJ 1/2 - (3-1) 

The bounds on the velocities v (of p. 310 and v a (of (|1.5p ) are a direct consequence of the Biot-Savart 
law and the conservation of the L 1 and L°° norms of vorticity. 

Proposition 3.2. Let q m (x) = qoXn in ( x ) be a multiple of the characteristic function of a simply connected 
bounded domain fl* n with C ,7 , 7 £ (0, 1), boundary. Then the global solutions v and v a of Euler and Euler-a 
equations (|1.3[) and (11.5[) . respectively, are in Lf oc (K; Lip (R 2 )) and for all a> 0, t £ R 

\\Vv(;t)\\ L ^<\\Vv m \\ Loo e c ^, (3.2) 
\\Vv a (;t)\\ L ~ < ||V« m || ioo e C|t| , 
where C = C (q m ) , independent of a, and the boundary of the vortex patch remains C 1 * 7 for all time. 
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For the Eulcr equations this result has been shown in [10]. We adopt their proof to show the uniform, in 
a, bound on ||Vi> a (-,i)|| L30 . The vortex-patch problem is reformulated in terms of a scalar function ip a (x, t) 
that defines the patch boundary by (t) = { x G R 2 |(^ Q (x, t) > 0} and is convected with the flow by 



dtp a 
~dt 



+ {u a • V) <fi a = 0, 



<p a {x,0)=<p in (x). 

To apply the method of the proof used in [10], the only ingredient we need is to show that Vu Q is uniformly, 
in a, continuous in the tangential direction of the boundary. First, we recall some properties of the kernel 

K a 



K a (x) = V ± ^ a {\x\) = y-rD^ a (\x\) , 

\X\ 



where 



b(£)+Iog. 



(3.3) 



(3.4) 



d^ a 1 

WW . _ W . - 



a \a/ r 



1 



D^ a (\x\) = - 
D H<* (r) = i- 



1 



a \x\ 



ar \a/ a z r \a/ a 6 \a/ r J 



The functions K and K\ denote the modified Bessel functions of the second kind of orders zero and one, 
respectively. For details on Bessel functions, see, e.g., [77]. Derivatives of decay to zero as ^ — > oo; and 
as — — > satisfy 



D* Q (r) = 



1 r , r ^ ( r \ 
log-+0 — , 

a \ a / 



4-7T cr 
1 1 



(3.5) 



D z m a (r) = -- 2 lo S- + ° -2 



47r cr 



1 



a 



2 ; ' 



D 3 * a (r) = -- — — + 



1 1 

4-7T ra 2 



log - J , 

a/ 



where the constants in the big O are independent of a. The filtered velocity gradient is given by 

Wu a (x,t)= f WK a (x-y)q a (y,t)dy 
Jr. 2 

As in [10], denote by W a divergence free vector field which is tangent to dQ, W = \7 ± tp. Then 
Wu a (x) W(x)= f VK a (x - y) q a (y) [W (x) - W (y)} dy, 

JR 2 

and we have the following result corresponding to Corollary 1 and Lemma in the appendix of [10]. 
Lemma 3.3. For 7 e (0, 1) 



\Vu a W\\ co „<C\\Vv a \\ LOO \\W\\ co „. 
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Proof. We stress that all the constants C are independent of a. We write 
Vu Q (x) W (x) - Vu a {x + h)W{x + h) 

^ VK a (x - y) q a (y) [W (x) - W (y)} dy 

VK a (x + h-y)q a (y) [W (x + h) - W (y)} dy 

VK a (x - y) q a (y) [W (x) - W (y)} dy 



\x-y\<2\h\ 



VK a (x + h-y)q a (y) [W (x + h) - W (y)} dy 



x—y\<2\h\ 



VK a (x - y) q a (y) [W (x) —W(x+ h)] dy 

x—y\~>2\h\ 

[VK a (x-y)- VK a (x + h- y)] q a (y) [W (x + h) - W (y)] dy 

x-y\>2\h\ 

h + h + h + h- 



Using the fact that \\7K a (x)\ < (cf. 



), we obtain that |2i| , \I 2 \ < C \\q a \\ L ^ ||W|| C0> -r h" 1 , also 



due to \D 2 K a (x)\ < -r^ (cf. (O, (El), we have |7 4 | < C (7) \\q a \\ L °o ||W|U ^ W . To bound the term J 3 

I I \X I 

we consider two cases \h\ < a and \h\ > a separately. We have 

\h\ < \\W\\co„ h-r\J\, 
where J = Si x _ y \ >2 \h\ Vif Q (x — y) q a (y) dy. First, let \h\ < a, write J as 

J= [ VK a (x-y)q a (y)dy- f V ' K a (x - y) q a (y) dy 

JR 2 J\x~y\<2\h\ 



Vu a (x) - Jl 



By dS3])-(IS31), we obtain 



\Ji\ < h a \\ L ™ I \VK a {x-y)\dy 

J\x-y\<2\h\ 



< \\Q a \\ L ~ 
<C\\q a \\ L 



<-<2 



1 1 

47r a 2 



log- 



y\ 



Now, let \h\ > a, then we write 



J = 



\x-y\>2\h\ 

= J2 + h- 



(VK a (x-y)- VK (x-y)) q a (y) dy - 



VK (x-y) q a (y)dy 



x-y\>2\h\ 



A bound for | J3I is obtained by using a lemma due to Cotlar (see [70] p. 291) 



Us I <c 



p.v. / VK (x - y) q a (y) dy 



C(\Wv a \\ L ^ + \\q a \\ Lx ) 
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We bound |J 2 | using the facts that if ,ifi > satisfy J R2 ±K a (^j dx = 2tt, J^ >2a (^j dx 

ttK (2), and we obtain 



\J2\<C\\q a \\ L 

<C\\q a \\ Lx . 
To conclude, in both cases we have 



\x— y\>2a 



1 



a\x-y\ 



\x-y\ 



a 



\k 



\x-y\ 



dy 



\h\<C\\W\\ COn h^(\\Vv a \\ L „ + \\ q a \\ LX ). 



□ 



Next we briefly recall the definition of homogeneous Besov spaces. Let S (W 1 ) be the Schwartz space and 
denote by Z' (M™) the dual space of Z (K") = |/eS (M«) : D 13 f (0) = for every multi-index f3 e N n |, it 

can also be identified as the quotient space of S /V, where V is the collection of all polynomials. Here (p and 
T '((/?) denote the Fourier transform of p in M™. 

We recall the Littlewood-Palcy decomposition. Choose a radial function ip e S such that <p € (R n \ {0}) 
satisfies supp^ c {§ < |£| < §} and J2^L-oo <Pq (£) = 1 f° r £ 7^ 0> where ip q (£) = </3(2~ 9 £), that is, 
<p, (x) = 29"(^(2«a;). For q G Z one defines the dyadic blocks by A q f = T^ 1 [T (<p q ) T (/)) = v? g * /. 
The formal decomposition / = ^2 q *L_ ao A<?/ n °ld s true modulo polynomials and is called the homogeneous 
Littlewood-Palcy decomposition. For set, 1 < r < oo, the homogeneous Besov space is defined as 



Sir = {/ G Z' (R n ) : 



< oo 



}■ 



where ||/||^ r = (E 9 °°=-co 2"' II AJ||^) ^ for r e [l,oo) and 
detailed description of Besov space see, e.g., [8,51,64,71]. 
We use the following results: 



Lemma 3.4. [60] For < (3 < 1 



liv/lli^. 



sup 9eZ (2 s "||A 9 /|| i2 ). For the 



(3.6) 



Lemma 3.5. Let seR, /3 > 0. T/iere exist constants c and C such that 



3 +3 < 



+0 



(3.7) 



Proof. This lemma follows directly from definition of the Besov spaces since 

A q ((-A)* /) || l2 = I^K) ^ ((-A)* /) || i2 = \\<p (2-^) |^ / (0 
and because (2"«£) is supported in {^e K" : §2« < |£| < §29}. 
Lemma 3.6. Let < (3 < 2, g e i?f and / = (l - c^A) -1 5 . Tften 



□ 



CT 



(-A)^ /2 / 



(3. 
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Proof. We have 



A q ((-Af /2 f 



L- 



L 2 



= a? 



L- 



(l + a 2 \i\ 



1 + a 2 \iX 

p ^ 1 2 

— (0 {2- q 0HOYdz 



a 2 l£l 2 



< 



sup ■ 



/3 



v y>o (l + yY 
< \\\g\\ L 2- 



1/2 



From which (13.81) follows 



□ 



The two following propositions show that the vortex patch vorticity, with C 1 ' 7 boundary, evolving under 



1/2 

Euler-a equations is in a homogeneous Besov space B. 



2,oo- 



Proposition 3.7. [60] If ft is a C ,7 , 7 > 0, bounded domain, then xn S B 
Proposition 3.8. Let v <E L°° ([0, T] ; Lip (R 2 )) mi/i divw = 0. Let tp solve 

^ + (v • V) p = 0, 
^(x,0) = ^"(x), 



3I/2 



with ip m e B]!^. Tin a 



,1/2 



< \\<P %n \Um e 



Cf£\\Vv(- 



- (It 



This result is a straightforward adaptation of the Proposition 3.1 of [1] to the homogeneous Besov spaces. 
Propositions 13.21 13.71 and 13.81 imply the following result: 

Corollary 3.9. Let q m (x) = qoXn in i x ) be a multiple of the characteristic function of a simply connected 
bounded domain Cl m with C 1,7 , 7 > 0, boundary. Then the global solutions q and q a of Euler and Euler-a 
equations (|1.3[) and (| 1 . 5j) , respectively, satisfy for all a > 0, t £ R + 



\\q(;t)\ 



< Ik™ II -i/2 e c -/'oll v ' u (-' 7 ")^-' ir , 

" "-^2,00 

< ||<f"|| -i/2 e^ollV'^-^Hi,-* 



Next we state a lemma that play an important role in estimating the convergence rate. 

Lemma 3.10. Let q' m — qaXn in > with Ct m being a simply connected bounded domain Cl m with C ,7 boundary, 
7 G (0, 1). Then for all T £ [0,oo) the solution u a of the Euler-a equations Ij 1.4)1 satisfies 



/2 



(3.9) 



for allte [0,T]. 



14 



Proof. Let < < \. Since (-B^f) = B% >00% for /3eR (see, e.g., [71]), we have 

c/ +1/2 \\Au a f L2 = a p+1/2 f Au a (x) • Aw a < Cc/ +1/2 ||Au Q ||^ ||Au a |U^ 



We stress that all the constants C are independent of a and /3. Now, by (|3.8[) and (|3.7[) 

^ +1 / 2 ||Au Q ||^ 



a^/ 2 || A ^|| < 



(-A) J 







(-A) 1/2 v a 




*f,oo " 





and by p. 61 



: ? <C||W||£ a \\Au°\\ l L 



1-/3 



< C\\q m \f L2 || A;/ 



a||l-/3 
L 2 ' 



Summing up we have 



hence 



lL 2 



lL 2 ) 



a /3+l/2 || Au «||» a < C \\q<*\\^ IIA^" 1 "^ 



||A^||, 2 <C||g«||^||^||2 (1+ «, 

2 , oo 

and taking the limit as (3 — ► we obtain 

a*||Au^<C||3%v ! i • 



□ 



In the next theorem we show that the solution of the Euler-a equations differs from the solution of the 
Euler equations by order (a 2 ) 3 ^ 4 , both having the same vortex patch as initial data. 

Theorem 3.11. Let q m = (joXn*™; where Xn in * s a characteristic function of a simply connected bounded 
domain Cl m with C 1 ' 7 , 7 E (0,1), boundary, and let v,v a be the solutions of the Euler and the Euler-a 
equations (|1.2[) and (II .4[) . respectively, with initial data q m . Then the difference w a = v — v a is square- 
integrable and obeys the estimate 

IKM)H 2 L2 <a^ 2 CV(t)\\q m \\^ /2 e cv ^\ 



where 



V(t)= (||Vz;(-,t)|| Loo + ||Vu q (. 



,)dr. 



In particular, there exists a constant C = C (q m ) such that 



\\w^;t)\\ L2 <a^ 2 Ce^ C \ 
Proof. The difference w a = v — v a satishes the following equation 

dW + (w a • V)« + VjVwf + {u a ■ V)w a + w^Vuf 



i)t 



(a 2 Au a ■ V)v - VjV (a 2 Au^) - VjVvj +V{p-p a ) = 0. 



(3.10) 



We remark that v and v a , obtained from the vortex patch vorticity by convolution with the Biot-Savart kernel 
are not in L 2 , however, the difference w a E L°° ([0, T] , L 2 (M 2 )) , since curlw Q is compactly supported and 



15 



/ K 2 curl w" (x, t) dx = L 2 (q m ~ q m < a J (x) dx = (this could be seen by using an asymptotic expansion of the 
kernel, see, e.g., [57], p. 321). Also, all the terms in (|3.10[) are in L°° ([0,T] , L 2 (R 2 )), since the velocities are 
in L°° ([0, T] , L°° (K 2 )) , and their gradients, as well as gradients of the pressures, are in L°° ([0, T] , L 2 (M 2 )) . 
We take an L 2 -inner product of (|3 . 1 0[) with w a (x, t) and obtain that 



2dt 

where 



onlK(-,*)III 2 <h + h + h, (3.H) 



h 



h = a 2 \{{Au a ■ V)v,w a )\ + a 2 \((w a • V)w,Aw Q )|, 
/ 3 = |(V (p- P a ),w a )\. 

here we used the identity, for h divergence free, 

((/ • V) ff + gjVfj,h) = ((/ • V). 9 , h) - ((h ■ V)g, f) . 

The first term is estimated by 

h<\\w a \\ 2 L2 \\Vu a \\ Lx . 

For the second term by (|3.9[) we obtain 

I 2 <a 2 \\Au a \\ L2 \\Vv\\ Lao \\w a \\ L2 

It remains to estimate the third term. We remark that the pressure is determined uniquely up to a constant, 
and in order to ensure that p £ L 2 , we require a side condition J R2 p (x) dx = 0. Notice that the mean free 
pressure of Euler equations satisfy 

/' = (viVj) 

where Ri — (— A) -1 ^ 2 ^ is the Riesz transform (see, e.g., [34]), hence by the properties of the Riesz 
transform and the Sobolev embedding theorem we have 

\\p\\ L , = \\WiWv < Ml* < C \\v\\ 2 wlA/3 <C\\q in f Ll/3 . 
For the Euler-a equation we have 

- <° ? <> + i a? <"? )2 - k (° 2a "?s-»?) < 

we write 

p a =p?+p%, 

where the mean free pressure pf , J R2 Pi {%) dx = 0, satisfies 



d d , „ . Id 2 



and pf satisfies 

As described above for the pressure of the Euler equations, we have 

,2 



M\\»<C <f 
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<a 2 ||Au«|| i2 ||Vw Q 



L- 



For the we have 

l|V^|| i2 = 
Hence, using also p.9p , we obtain 

/ 3 <||V^|| i2 \\w a \\ L , <Ca 3 / 2 \\q a \ 

From the above and p. lip , and by using Gronwall lemma, we obtain 

IK (;t)\\ L2 < Co?l 2 \\q m \\ m V(t)e CV{t \ 

whereF(i) = J* (|| Vw (•, r) \\ Loo + || Vu a (•, t) || Loo ) dr. Hence by ([321) 

IKM)|| L2 <Ca 3 / 2 ||g m ||^ /2 ||V U m || ioo e c 'e C l! V ^ll 
<a 3 / 2 Ce CeCt , 
where C depends only on the initial data q m . 



□ 



Appendix 

A.l Global regularity of Contour Dynamics-a equation 

We consider the vortex patch problem, i.e., a system in which the initial vorticity q m is proportional to the 
characteristic function of a bounded domain Q, m , q m — qoXfi in i under the evolution of the Euler-a equations 
(|1.5p . Due to the conservation of the two-dimensional Euler-a vorticity along particle trajectories for such an 
initial data (in fact, it is enough for the initial vorticity to be in the space of Radon measures in R 2 , see [63]), 
the vorticity q a (t) remains a characteristic function of an evolving in time domain Q a (t) . In this section 
we present the result which states that the boundary of a vortex patch evolving under Euler-a equations 
(jl.5p remains as smooth, for all time, as initially boundary, provided the latter is smooth enough in a sense 
specified in Theorem I A. 1 1 below. 

In two dimensions the evolution of the boundary x (a, t) of a vortex patch under the Euler-a equation 
(|1.5p is given by 



dx , . 



x (a, 0) = x m (a) , 



I * a (\x(a,t)-x(a',t)\)^(a',t)da', 
l s i da 



(A.12) 



where 



* a 0) = — 



Koi- 



logr 



and Ko denote the modified Bessel functions of the second kind of order zero. As we mentioned before, 
for details on Bessel functions, see, e.g., [77]. The integro-differential equation (|A. 12|) is an analogue of the 
so-called contour dynamics (CD) equation [57,79], describing the evolution of the boundary of vortex patch 
evolving by means of the Euler equations (JTT 



(A.13) 



Tt (CT ' l) = ~fr J 1 l0g {]X (<T ' t] ~ X (<T '' m t} df7 '' 

x (a, 0) = x m (a) . 

The CD-a equation (|A. 12|) is derived following arguments similar to those used in the Euler case, see [57,79 
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We show that CD-a equation (|A.12[) is well-posed in the space of Lipschitz functions and in the Holder 
space C n, @ , n > 1, which is the space of n-timcs diffcrcntiablc functions with Holder continuous n th derivative. 
Let us first describe the Holder space C"'' 3 (S 1 C R;R 2 ),/3 £ (0, 1], which is the space of functions x : S 1 C 
R — > R 2 , with a finite norm 

J\ d k d" 

where 



n,/3 - 

k=Q 



d k 



CCS 1 ) 

IMIc"(si) = SU P \ X ( a )\ 



da nX 



and |-L is the Holder semi- norm 



o-es 1 

\x(a)-x(a')\ 
a^'es 1 \a-a'f 

The Lipschitz space Lip (S 1 ^ is the C 0,1 space, that is, with the finite norm ||^|lLi p (si) = H^llc^s 1 ) + \ x \v 
We also use the notation 

_ \x{a)-x{a')\ 

We consider the CD-a equation (|A.12[) as an evolution functional equation on either one of the following 
Banach spaces: Lip, C 1 -' 3 , (3 S [0, 1], C 2 ' 7 , 7 S [0, 1), or C n ' e , n > 3, 6» e (0, 1). We have the following result 

Theorem A.l. Let V be either one of the following spaces: Lip (S 1 ) , or C 1 ' 13 (S 1 ), (3 G [0, 1], or C 2 < 7 (S 1 ), 
7 £ [0, 1), or C n ' e n > 3, 9 e (0, 1). Let x m £ V n > 0}, then there exists a unique solution 

x £ C 1 ((— 00, 00) ; V fl {l^l^ > 0}) 0/ (|A.12[) rofi initial value x(a,0) = x m (a). Ln particular, if Xq £ 
C°° (S 1 ) n {|a:|, > 0} then x £ C 1 ((-ex, 00) ; C°° (S 1 ) n > 0}) . 

We remark that, although the kernel ^> a and its first derivative D^ a are continuous bounded functions, 
its higher derivatives D nll fy a 1 m > 2, are unbounded near the origin, and the chord arc condition |a;|* > 0, 
which implies simple curves, allows us to show the integrability of the relevant terms. 

We only sketch the main steps of the proof, since it is in the spirit of [4,5] and [57, Chapter 8], which can 
be consulted for details of such a proof. In [5] we show the well-posedness of vortex sheet problem for Euler-a 
equations, and it contains various estimates involving the derivatives of the kernel ^f a , and [57, Chapter 8] 
describes the proof of the original vortex patch problem in the Euler equations case. 

The following are the main steps involved in the proof of Theorem IA.1I In the first step, we apply 
the Contraction Mapping Principle to the CD-a equation (j A. 12|) to prove the short time existence and 
uniqueness of solutions in the appropriate space of functions. Next, we derive an a priori bound for the 
controlling quantity for continuing the solution for all time. At step three one can extend the result for 
higher derivatives, using the estimates derived in [5] and [57, Chapter 8]. 

A. 1.1 Local existence of Contour Dynamics-a equation 

Next we show the local existence and uniqueness of solutions in the Lipschitz space, the details and the 
estimates in other appropriate spaces can be done in the same spirit following arguments presented in [5] 
and [57, Chapter 8]. First we recall some properties of the kernel $ a , see also ()3.4[) - (|3.5[) . For ^ — > 

* Q W = log" + 0(1), (A-14) 

and for large — 

* a (r) =OQogr). (A. 15) 
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Also, D^ a is bounded for all r G [0, oo), 



DV a (r) = O [ - ) , 



(A.16) 



where in the big O the constants are independent of a. To apply the Contraction Mapping Principle to the 
CD-a equation (|A.12[) we first prove the following result: 



Proposition A. 2. Let 1 < M < oo, and let K be the set 



K M = jz G Lip (S 1 ) 



x\\u P < M >\ x \. > 



M 



Then the mapping 



x(T)^u (x (T))= / K a (x (r) - x (r') ) dV 
Js 1 

'M 



(A.17) 



defines a locally Lipschitz continuous map from K , equipped with the topology induced by the |H| Lip norm, 
into Lip. 



Proof. We start by showing that u {x (a)) maps K M into Lip. Let x G K M . By fA~T4|l and (fA~T5)) we have 



\u(x(a))\<q / * a (\x(a)-x(a')\) 
Js 1 



da 



da' 



(A.18) 



^CgoCloga + l + logM)!^. 



To show Lipschitz continuity of u (x (a)) we use that by mean value theorem and (I A. 16|) . we have that for 
x (a") G B (x (a) , \x (a) — x (a)\), the ball centered at x (a) with the radius \x (a) — x (<r)|, 



|* Q (\x (a) - x (cx')|) - * a (\x (a) - x (a')\)\ < {\x (a") - x (a')\) \x {a) - x (a)\ 

C 

< — \x(a) -x{a)\ , 
a 



(A.19) 



and hence 



|u (x (a)) -u(x (a))\ < q Q / \* a (\x (a) - x (a')\) - # Q (\x (a) - x (a') 

Js 1 

C 

< —qo Mi \x{a) - x(a)\ . 



da 



da' 



Now, we show that u (x) is locally Lipschitz continuous on K . It is enough to prove that for x G K , 

y e Lip (s 1 ) 

(A.20) 



\\D x n(x)y\\ Up <C[-,M,\\x\\ 



Lip I liyilLip • 



Let x G K M , y G H 1 (S 1 ) , we compute 



D x u (x {a)) y (a) = —u (x (a) + ey (a)) 



e=0 



- qo f mr (\ x (a) - , (.01) {x {a) ~ * ( f? • v (g0) £ V) * 

J S i \x(a)~x{a')\ da 

q f *°(\ x (a)-x(a')\)^(a')da' 
J s i da 

:F 1 (x(a))y(a)+F 2 (x(a))y(a). 



19 



Next we show (|A.20|) . To estimate the C° norm we use (|A.16[) for F± 



|fiW<r))l/WI<« / DV*(\x(a)-x(a')\)\y(a)-y(*')\ 
Js 1 



do-' 



C 

< —9o ||y|| c o ■ 



and (Q4|> and (jA~l~5l) for F 2 



\F 2 (x(a))y(a)\<q / |* Q (\x (a) - x (a') 
Js 1 



dy , „ 

Ta^' 



da' 



<Cq (loga + l + logM)^. 
Next we show Lipschitz continuity of D x u (x (a)) y (a). For F2 one uses (|A. 19|1 . For F\ we have 
\F 1 (x(a))y(a)-F 1 (x(a))y(a)\< 



< q / m a (\x (a) - x (o-')l) - D*" (\x (a) - x (o-')|)| \y (a) - y {a')\ 
Js 1 



da 



da' 



+ 2,o Z?* (|x (a) - x (a )|) |x (a) - * (a)| 



da' 



+ qo / D* a (\x(*)-x(a')\)\y(a)-y(a)\ 
Js 1 

= h+h + h- 



dx 
da 



(<0 



da' 



For Ii, by the mean value theorem, (|3.5p and due to the fact that > jj, we have that for a" S S 1 and 
such that x (cr") E B (x (a) , x (cr) — x (a)\) 

\D^ a (\x {a) - x (a')\) - D^ a (\x (a) - x (a')\)\ < \D 2 V a (\x (a") - x (a')\)\ \x (a) - x (a)\ 



11 

Air a 2 



log 



\x(a")-x(a')\ 



1 



< C{M)^\a -a 



log 



a — a 



h < \a — a\C (M) —qo \\y\\ c o \x\ 1 



Therefore, 



<C{M,-,q ) \\y\\ c0 \a-a\. 
For I 2 and I 3 we use (|A.16[) and \x\^ > jt to obtain 



log 



1 da' 



h,h<C (M,-,q )\y\ 1 \a-a 



a 



□ 



Proposition I A. 21 implies the local existence and uniqueness of solutions: 

Proposition A. 3. Let K M = jx e Lip (S 1 ) : ||x|| Lip < M, |x|, > jj^and let x G Lip (S* 1 ) n {\x\„ > 0}, 

then for any M, 1 < M < 00, such that xo E K , there exists a time T(M), such that the system (|A. 12|) 
has a unique local solution x E C 1 ((—T(M), T(M)); K M ). 
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A. 1.2 Global existence of Contour Dynamics-a equation 



To show the global existence of the CD-a equation, we assume by contradiction, that T max < oo, where 
[0, T max ) is the maximal interval of existence, and hence the solution leaves in a finite time the open set 



K , for all M > 1, that is, limsup^y- 



"llv = 00 or limsu Pt^T- ax H" -,t)|, 
global bounds on ra^ji and ||x (•, t)\\ v in [0, Tmax), we obtain a contradiction to the blow-up, and thus the 
obtained local solutions can be continued for all time. The result extends to negative times as well. 



= oo. Therefore, if we show 



To control the quantities 



and ||a;(-,i)||y one needs to bound J Q 



\x(-,t)\, ii* V' "JHV 
proposition shows the bound on ||Vu Q || L=0 for the vortex patch initial data. 



\S7 x u a (x(-,t),t)\\ La 



dt. Next 



Proposition A.4. Let q m e L 1 (R 2 ) n L°° (l 



Proof. We write 



Using (|3. 5j) we obtain 



and 



\\V x U a (X(;t),t)\\ 



then 



< C 



(Ik* 



IL 1 / 



V x u a (x,t)= \7K a (x-y)q a (y,t)dy 



R 2 

+ / =h+h- 

\x — y\<a J \x— y|>ct 



h<C\\q a (;t)\\ Lao r\D^ a (r) 
Jo 

<C\\q m \\ L; 



rdr 

1 1 , r C 
log - + — 

47r a a z 



dr 



< c\W 



h < sup \VK a (x-y)\ f 

\x-y\>a JR 5 



q a (y,t)dy 



< C- 



IL 1 



□ 
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